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ABSTRACT 

The sample of 137 low-redshift type la supernovae with 0.05 < z < 0.3 obtained 
from the SDSS-II Supernova Survey for the southern equatorial stripe of 300 square 
degrees is used to derive the luminosity functions of type la supernovae and of their 
host galaxies in the gri passbands. We show that the luminosity function of type la 
supernova host galaxies matches well with that of galaxies in the general field, suggesting 
that the occurrence of type la supernovae does not favour a particular type of galaxies 
but is predominantly proportional to the luminosity of galaxies. The evidence is weak 
that the supernovae rate varies with colour of host galaxies. The only evidence that 
points to possible correlation between the supernova rate and star formation activity is 
that the supernova rate in late-type galaxies is higher than that in early-type galaxies 
by 31 lb 35%. In our low redshift sample the component of type la supernova rate that 
is proportional to star formation activity is not manifest in the integrated supernova 
rate, while our observation is compatible with the current two-component models. The 
sample contains 8 type la supernovae whose host galaxies were not identified, but it is 
shown that their occurrence is consistent with them occurred in low luminous galaxies 
beyond the survey. The luminosity function of type la supernovae is approximately 
Gaussian with the full-width half maximum being a factor of 1.4 in luminosity. The 
Gaussian distribution becomes tighter if the ratio of extinction to reddening, Ry, is 
lower than the characteristic value for the Milky Way and if luminosity is corrected for 
the light curve shape. The colour excess is ~ 0.07 mag which is significantly smaller than 
reddening expected for field galaxies. This colour excess does not vary with the distance 
of the supernovae from the centre of the host galaxy to 15 kpc. This suggests that the 
major part of the colour excess appears to be either intrinsic or reddening that arises 
in the immediate environment of supernova, rather than interstellar reddening in host 
galaxies, and most of type la supernovae take place in relatively dust free environment. 



Subject headings: supernovae : general 
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Introduction 



Recent studies have revealed that our understanding of the mechanism of Type la supernovae 
(SNe la) is poorer than we had thought. It is suspected that there are two different progenitor types 



(Dallaporta 


197.-^: 


Tammann 


1982: 



20061 ). The evidence, however, is not conclusive yet, and some observations at high redshift do not 
fit this picture; for example, there are indications for a drop in the S N rate at z > 1.5, where star 
formation rate is still rising (jPoznanski et al.ll2007l : iDahlen et al.ll2008l ). which does not support the 
presence of a large prompt component. 

The present paper studies the luminosity functions (LF) of S Ne la and their host g alax- 



ies and correlations in properties of SN la and of host galaxies. iMiller &: BranchI (jl990l ) and 



Richardson et al 



(12002) stu died the LF of various types of SNe bas ed on the Asiago Supernova 



Catalog (IBarbon et al.lll989l ). In particular, iRichardson et al.l (120021 ) showed that the LF of SNe 
la is consistent with the Gaussian distribution with (Mb) = —19.46 + 51og(i?o/60) and a = 0.56 
using 111 spectroscopically normal SNe la without correcting for the light curv e shape parameter. 
With the correction for the decline rate a can be reduced to 0.11 (|Phillipslll999l ). We are not aware 
of a study of the LF of SNe la host galaxies. We may hope that a comparison with the LF of field 
galaxies may hint us as to what type of galaxies would preferentially host SNe la. 

We have accun iulated a sample of SNe la acquired in the second phase of the Sloan Digital 
Sky Survey (SPSS: lYork et al.ll200d ). Supernovae were searched during September to November 
of 2005 to 2007 by repeated imaging for the sky area of 300 square degree o f the southern equa- 
torial region, -60° < R.A. < +60°, -1?25 < Decl. < +1?25, every two days (jPrieman et ahlboosl : 
Sako et al.l boosl : IPildav et allboosl l. Approximately 500 SNe la at 0.05 < z < 0.4 were iden 



tified using w ell-defined selection criteria and their light curves were measured (jSako et al.ll2008l : 
Holtzman et a l. 200 8.) • The advan t age of the use of th e SDSS is, apart from its accurate five colour 
photometry (IFukugita et al.l 119961 : ISmith et al.ll2002l ). that most of the galaxies that host super- 
novae have already been photometrically observed with homogeneous presetted criteria, so that 
one can study the properties of those galaxies and correlations between the supernovae and their 
host galaxies. The disadvantage is that the SDSS is somewhat too shallow for this purpose due 
to a limit ation arising frorn the t ime-delay-and-integrate mode imaging with the 2.5 m aperture 



telescope (iGunn et al. 



1998 



20061 ). 



The SN candidates are spectroscopically fo llowed-up with oth er telescopes as much as the time 
allows to determine their types and redshifts (jZheng et al.ll2008l ). After the SN light faded away 
spectroscopic observations are carried out for its host galaxy to detremine redshift of the SN la 
candidate. In this paper we use the data from the first year for the completeness is high for redshift 
of host galaxies for the first year SNe. We limit the sample effectively to z < 0.3 for the sample 
completeness. Five colour photometry had been made for host galaxies. A O.IL* galaxy at z ~ 0.3 
would give g ^ 22.7, r ~ 21.9, i ~ 21.4, z ~ 21.1, so that we can sustain reasonable accuracy in 
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SDSS photometry for the four passbands (jHogg et al.ll200ll : llvezic et al.ll2004l ). For the u passband, 
however, a O.IL* galaxy at z ~ 0.3 gives u ~ 24.3, which is significantly beyond our survey limit 
u = 22.0, a calibre of our survey with a 2.5-m telescope, and hence, the u passband information 
cannot efficiently be exploited. 



In this paper, we take Hq 
parameters. 



70 km s ^ Mpc ^, = 0.3 and 17 a = 0.7 for the cosmological 



2. Sample 



In the 2005 SDSS-II observation, 706 SNe la candidates were identified including objects with 
observations at only a few epochs and 130 SNe la among them are spectroscopically confirmed 
with additional 16 identified as spectroscopically probable SNe la. The latter are the objects whose 
spectrum do not show significant Si II absorption-line feature while the overall spectra match with 
the la templet. Follow-up observations were not complete due to the time limitation while SNe 
candidates would be observable. We are left with the candidates whose light curves indicate SNe 
la but were not spectroscopically confirmed. Our experience with the spectroscopically confirmed 
sample tells us that spectroscopic fai lure rate to confirm SNe la is only 10% of the light-curve 
selected candidates (jSako et al.l l2008l ). This significant failure rate, however, applies to loosely 
selected candidates raised for spectroscopic follow-up that use light curves only in early days from 
the first detection, and it is substantially reduced for tightly light-curve selected SN la candidates, 
such as those used in this paper. A contamination of a few misidentified SNe la out of a hundred, 
if any, would not spoil global statistical analyses, such as those presented in this paper, beyond 
our statistical errors. Host galaxies of those candidates were spectroscopically observed later and 
their redshifts were measured. We refer to these cases as 'host-z measured SNe la'. There are 168 
objects in this category, which make our total sample to comprise 314 SNe la with known redshift. 

For the supernova sample, g, r, and i-band light curves are measured by fitting images with a 
model of the galaxy background and the SN point spread func tion without pixel re-sampling frame 
by frame (scene-modeling photometry) (iHoltzman et al.ll2008l). These SN m ulti-band photometric 
data are fit simultaneously by the SALT2 light curve fitter (iGuy et al.ll2007l ) to derive photometric 
parameters. Some SNe la in our sample are not suitable for detailed analysis: for SNe la discovered 
in the early or late phase of the observing season, i.e., early September and late November, light 
curves are observed only partially, with which it is often difficult to determine the accurate light 
curve shape and hence peak brightness. Some SNe la are faint and their signal-to-noise ratio at 
pe ak brightness i s too low due to their large distances. We set the following criteria, as were done 
bv iDildav et al.l (|2008l ). for the secure acquisition of light curves: (i) at least one measurement 
with Trest < —2 days; (ii) at least one measurement with Trest > +10 days; (iii) at least five 
measurements with —20 < Trest < +60 days; (iv) at least one measurement with signal-to-noise 
ratio above 5 for each of g, r, and i; (v) reasonable light curve fit whose reduced calculated by 
SALT2 is less than 3 for probable and host-z SNe. For spectroscopically confirmed SNe la, the 
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criteria was not applied. Here Trest is the rest-frame time from the epoch of peak brightness in 
the rest-frame -B-passband. 222 SNe la among 314 passed these li ght curve criteri a. The known 



pecu liar SNe la, SN2005gj (|Aldering et al.1 12009 : IPrieto et al.ll2007l l and SN2005hk (|PhilliDs et al 
20071 ) are rejected here. We limit the sky area to —50.0 < a < +55.0 and —1.25 < 5 < +1.25 and 



the time coverage to 53626 (13 September 2005) < tmax(MJD) < 53691 (11 November 2005) to 
avoid the edge effects. This leaves 207 SNe la with us, the distribution of which is shown in Figure 
[Has dotted histogram. We take this as the basic SN la sample. These selections, up to item (v), 
are readily implemented in the analysis procedure. The effective area of the survey is 262.5 square 
degree. 

To make the sample incompleteness well defined for distant faint supernovae and the sample 
suitable to an application of the 1/T^ax method to compute the LF, we make the sample magnitude 
limited. Figure [2] shows apparent maximum brightness in the r-passband as a function of redshift, 
which suggests us to set the limiting magnitude to be = 21.5mag. The figure shows that 
this ensures reasonable completeness of SNe la at z = 0.20 — 0.25. After applying this limiting 
magnitude, 137 SNe la are left with us for a magnitude limited sample, with the mean redshift {z) = 
0.207, which we take as our final sample for our analysis. Among magnitude limited sample, 72 are 
spectroscopically confirmed, 5 are spectroscopically probable and 60 are host-z SNe. Summary of 
the sample selection is shown in Tabled) 

The solid histogram in Figure [T] represents the redshift distribution of this magnitude lim- 
ited sample. The shade indicates the distribution expected for the empirical redshift distribution 
calculated using the code to generate llight curves SNANA0 (R. Kessler et al. 2009, private com- 
munication, submitted to PASP) assuming the SN la r ate, rsNeia. {l + z Y'^^^'^ , which is inferred 
from the observed SDSS SNe la rate at low redshift ([Pildav et al.ll2008l ). The conditions of the 
SDSS observations, the software search efficiency and the light curve criteria are taken into account, 
and the normalisation is determined from the number of SNe below z = 0.15, for which redshift 
range SDSS SN observation is complete. The detail of simulation is described in Appendix. The 
shade stands for ±20% uncertainties of the normalization, corresponding to the Poisson error of 
the number of SNe la below z = 0.15. The SN frequency expected for no evolution of SN rate is 
shown by the dashed curve, which lies within the error indicated by shades. This figure shows that 
our sample is likely complete io z < 0.20. 

The sample incompleteness at higher redshift is caused mainly by spectroscopic targetting and 
is given by fitting the ratio of observed number of SNe in our basic SN la sample to the simulated 
number as seen in Figure [3) 



e{z) 



1.0 Z<Zc 

1.0 — {z — Zc) /Az Zc < Z < Zc + Az 
0.0 z> Zc + Az, 



(1) 



^http : //sdssdp47 . f nal . gov/sdsssn/SNANA-PUBLIC/ 
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where the best fit values are = 0.162 ± 0.029 and Az = 0.279 ± 0.041. 



For each SN la the nearest primary object that resides within 10 ar csec from the SN is identi- 



fied a s its hos t galaxy in Catalog Arc hive Server of SDSS Data Release 6 (jAdelman-McCarthv et al 



20081 : see also IStoughton et al.ll2002l ). with the identification visually confirmed. When misidentifi- 
cation is suspected, it arises mostly from deblended galaxies. Figure H] shows the apparent r-band 
Petrosian magnitude distribution of host galaxies for 207 SNe la before the magnitude cutoff with 
the dotted histogram. The solid histogram is the host galaxies for our magnitude limit sample. 
There are 9 SNe whose host galaxies are not identified in the basic sample of 207 SNe (8 SNe in 
the 137 magnitude limit sample) within ~ 10 arcsec from a SN. These hostless SNe la are listed in 
Table [21 which gives the upper limit on absolute brightness of host galaxies searched in our survey 
to be 0.025 — 0.09L* for the apparent magnitude limit of r = 22.2 mag. The physical distance 
corresponding to 10" is also given in the table. Petrosian magnitudes of u, g, r, i, and z-passband 
are use d to compute rest-frame absolute brightness of host galaxies in 5 passbands using kcorrect 
v4_l_4 (jBlanton &: RoweisI 120071 ) . Redshifts are fixed to the spectroscopic values. 



Application of the 1/Fmax method 



The light curves of SNe are fit by the SALT2 (|Guv et ahlboOTl ^. which yields apparent 



maximum 



brightness in the rest-frame i?-passband, niB, together with the parameter describing the spectral 
variation xi and the optical depth representing the colour excess at maximum brightness c. The 
parameter xi is the weight of the next leading component of the spectral templet and is related 
to the light curve shape parameter Amis (Phillips 1993), as Amis = 1-09 — O.lGlxi -|- 0.013xf -|- 



Oixf) (Guv et al. 20071 ). SALT2 corrects for Galactic extinction following the extinction map of 
Schlegel et al.l (|l998l ). 



We apply the l/Vmax method to calculate LFs of SNe la (and of their host galaxies) to the 
magnitude limit sample of 137 SNe whose maximum brightness is brighter than 21.5 mag in the 
r-passband in the observed frame, by calculating the maximum redshift, Zmaxi at which each SNe la 
(or their host galaxies) can be observed within the magnitude limit. Since the multicolour passband 
light curve fit by SALT2 gives us the spectral energy distribution of SN la, characterised with the 
shape parameter and the colour excess parameter, we can calculate expected apparent brightness 
in the r-passband in the observed frame at various redshifts using this spectral energy distribution. 
The limiting redshift gives the maximum volume that can be surveyed, 14nax! as 

where uj is the solid angle 105 x 2.5 = 262.50°, and the index i refers to SNe la. Since ^^x is 
regarded as the volume surveyed to find i-th SNe, the LF is obtained by summing the inverse of l^ax 
within specified magnitude bins (AM) assuming that the LF is not evolving over the respective 
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redshift range. Taking into account the effect of visibility time and spectroscopic incompleteness, 
the LF is calculated as 

<^(M)AM(Mpc-3) = ^ ( iy (3) 

ie\M^-M\<AM/2 ^^"^ ^ ' 

where ct* is the time of visibility in the rest-frame at which each SN la is observed. If the SN 
would be observed just at one epoch, as was done in a number of observations to derive the 
SN rate, the visibility time will be a time span over which each SNe can be detected above the 
detection limit: in this case a fainter SNe would have a shorter visibility time. In our case, however, 
observations have been made for the same field of sky continuously with the magnitude limit set 
for peak brightness. The visibility time will then be a time span of the survey observation. From 
the criteria on the date of maximum brightness, the visibility time in the observed frame is 65 days, 
and hence ci* = 65/(1 + z*) days in the rest-frame. For SNe la M* is absolute peak i?-passband 
magnitude of SNe la, whose apparent rest-frame magnitude is estimated from SALT2. For host 
galaxies it is absolute magnitude in the rest-frame estimated using kcorrect. The factor ei^z) is 
the completeness correction, eq. ([TJ. The LF is represented in units of per Mpc^, and r (year) 
absorbs the time of the duration of observability when we deal with SNe. We take r = 1 yr as the 
unit. We describe in Appendix simulations we made to show that sample incompleteness and our 
corrections do not induce particular systematic errors to our analysis and our procedures allow us 
to recover the true LF, SN rate and related quantities. 



4. Luminosity function of SNe la 



To estimate the intrinsic brightness of SNe la, we must correct for dust extinction within host 
galaxies. In fact, SNe la show the variation in colour that could be attributed to dust extinction 
within host galaxies and/or interpreted as an intrinsic colour variation. The colour information is 
obtained from the colour excess parameter c of SALT2, which is defined by c = [B — l^)max — ((-B — 
^)max) at -B-passband maximum brightness, where the second term is colour of the SN la templet. 
Figure [5] shows the distribution of this colour excess c = E{B — y) . The distribution is asymmetric 
with respect to E{B — V) = 0, similar to that of the SNLS sample (jAstier et al.ll2006l ). If we assume 
the color distribution as an exponential distribution oc exp(—E{B — V)/A) smeared by intrinsic 
Gaussian color distribution with the dispersion of a, the observed distribution can be fitted with 
A = 0.048 and a = 0.074 as in Figure El Th is value of A is smaller than A = 0.138 obtained for 
the nearby SN la sample by I Jha et al.l (|2007l ). T he mean a i id dis persion (c) = 0.176 ± 0.280 are 
compared with {E{B - V)) = 0.128 ± 0.173 from I Jha et aD (|2007l ). Several data points with high 
value of E(B — V) greater than 0.5 contribute significantly to c. If they would be removed the 
mean and dispersion will become (c) = 0.061 =b 0.107. One may suspect that our sample may lack 
highly extinct or very red SNe due to selection effects. It is unlikely, however, that we dropped 
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such SNe at least at lower z where SNe are surveyed deep enough to reach brightness significantly 
fainter than are dimmed by extinction, unless E(B — y) > 1 which are dropped as our selection 
procedure for the SN la sample. The fact that the distribution does not change appreciably from 
low z to our redshift limit indicates that the decline of the event rate that give E{B — > 0.2 is 
not due to the selection effect and SNe that give a larger E{B — is missing. The distribution of 
E{B — y) corresponds, if interpreted as extinction, to extinction in the V passband ranging from 
—0.14 mag to 0.53 mag for one sigma, if Ry = 3.1 is assumed. The range of positive numbers is 
consistent with what is expected from extinct ion in the galactic discs with the HI column density 
< 9 X 10^0 cm-2 (jBurstein h Heileslll978l ). which is a typical value for disc galaxies unless they 
are edge-on. Extinction may also receive a contribution from SN itself and its environment. We 
come back to this problem in the end of the next section. 

Figure [6] shows the distributions of Ay for the SNe la sample, taking the colour variation as 
reddening and assuming Ry = 3.1. Note th at above Ay = 3 SNe la candidates, if any, are rejected 



during SN identification with light curves (jSako et al.ll2008l ). This is compared with Ay for star 
forming field galaxies (dotted histogram) , as calculated from the Ha / H/3 Balmer line ratio ob tained 
by SDSS within 3 arcsec of the central region of galaxies for z < 0.2 (jNakamura et al.l 12004 ). The 
mean value of extinction of SNe la, {Ay) = 0.20 mag, is significantly smaller than that for star 
forming field galaxies, (Ay) = 1.10 mag. This suggests that most SNe la take place in regions 
clearer in dust than actively star forming regions. 

We provisionally assume that the colour variation is due to extinction within host galaxies and 
that the extinction law is the same as that in our Galaxy with Ry = 3.1. We then estimate bright- 
ness of SNe la by correcting for extinction as AIb — /3 x c with /? = Ry + 1 = 4.1. We note, however, 
that this ratio of total to selective extinction is n ot well .justified . Ther e are some indications that 
Ry for SNe la is lower than the Gala ctic value: lAltavilla et al.l (j2004l ) used the bluest SNe la to 



estimate intrinsic colour of SN la and Reindl et al 



(|2005l ) used SNe la i n early- type galax i es and 



far outlying SNe la in spiral galaxies. Both authors obtained Ry ~ 2.5. iNobili &: Goobaii (|2008l ) 
also found a low value of Ry = 1.75 simu ltaneously deriving templet s of SNe colo u r evo lution 



in time and extinction. lAstier et al.l (l2006l ) resulted in Ry = 0.57 and iKessler et al.l (|2009l ) gave 
Ry = 2.18. All Ry thus derived are smaller than the canonical value for the Milky Way. 

Figure [7] shows the LF of SNe la in the -B-passband, which is fitted with the Gaussian dis- 
tribution with the mean and dispersion given in the upper left corner of each panel. The panel 
(a) assumes that the colour variation c arises from extinction with /? = 4.1, giving the mean 

= —19.42 (in the Johnson zero point) and the dispersion 0.24 mag. If we would adopt a 
smaller value for Ry, the LF of SNe la becomes closer to a more regular Gaussian distribution with 
a smaller dispersion, as shown in Figure [T^b). The minimisation of dispersion with respect to /? 
results in /? = 2.93, which leads to the narrower Gaussian width of 0.16 mag with = -19.32, 
which is 0.10 mag fainter than with /3 = 4.1. Thus the smaller Ry gives more homogeneous SN la 
luminosity. 
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This small value of Ry obtained by minimizing the wi dth of Gauss i an fitt ed to SNe LF is an 
alternative manifestation of the low value of f5 obtained by lAstier et alj (120061) by minimis i nK the 
resi dual scatter in the H ubble diagram along with cosmological parameters. iKowalski et al.1 (|2008l ) 
and iKessler et al.1 (|2009l ) indicated that minimizing the scatter in the Hubble diagram tends to give 
Ry biased towards a value lower than the true value. Our simulation, however, shows that the 
value of /3 may be biased to the lower value but no more than by ~ 0.1, so that this cannot be the 
reason for Ry being significantly smaller than 3.1. We do not conclude here that Ry is actually 
smaller but take /3 = 4.1 as our fiducial choice, keeping in mind the uncertainty from Ry in the 
analysis in what follows. 

It has been argued that the maximum luminosity correlates with the light curve shape, or more 
specifically the decline rate (Amis, stretch, or SALT2's xi ), and the inclusion of the correlatio n 



with it makes the behaviour of the LF tighter (jPskovskiilll984l : iPhillips 



1993 



Hamuv et al.lll99ffl . 



We show in panels (c) and (d) of Figure [7] the LF where brightness is corrected by axi, with a 
chosen to minimise the width of the Gaussian. In panel (c), f3 is fixed to our fiducial value of 4.1 
and the minimisation gives a = 0.052. In (d) both a and f3 are chosen to minimise the width, 
which results in /3 = 2.52 and a = 0.123. The correlation of xi with brightness makes the Gaussian 
distribution narrower, especially for the case in which both a and /3 are optimised. The narrowest 
Gaussian is obtained with a = 0.123 which corresponds to a' = 0.76 where the correction for the 
light curve shape to brightness is expressed as a' Amis. This value is consistent with 0.78±0.18 
obtained by iHamuv et al.l (|l996l ) . 



Figure [8] shows the correlation between brightness of SNe with the shape parameter and the 
colour excess. This correlation is the reason why the width of the luminosity function decreases 
upon inclusions of the light curve shape parameter and the colour excess parameter (with a small 
Ry). The correlation in the upper panel is represented by Mb = —19.34 + 2. 93c as shown in Figure 
[TJb), and that in the lower panel by Mb — 4.1c = —19.43 — 0.052xi as in Figure [Tl^c). 



5. Host galaxies 

The LFs of SNe la host galaxies, as calculated by the 1/V^max method, are shown in Figures [H 
[To] and [TT] with solid histograms for the r, g, and i-passbands. Hostless SNe are indicated by shaded 
histogram at r i ghtmo st bin. We draw with solid curves the LF of general field galaxies obtained by 



Blanton et al.l (|200ll ). multiplied with the luminosity. The curves show good match of the LF of 
SN host galaxies with that of the field galaxies, which means that the LF of galaxies derived from 
SNe la faithfully represents that of galaxies in the field. We do not see any particular deviations 
between the two for three colour passbands g, r and i we studied, meaning that the occurrence of 
SNe la is primarily proportional to the luminosity of galaxies. Matching the two LF's gives the 
SNe la rate in the conventional supernova unit (SNu), the SN rate per lO^'^ solar luminosity per 
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century r^. Taking Mr{Q) = 4.62, we obtain 

TL = 0.227 ±0.027 SNu(r), (4) 



where we use the luminosity normahsation of iBlanton et al.l (|2003l ) shifted to z = 0.2 using their 



evolution prescription with the Q parameter that represents the evolution of luminosity per redshift 
interval^. This matching was done using data points at magnitude bins where the number of con- 
tributing galaxies are greater than 10 and error is Poisson from the SN number used in the l/l^nax 
analysis. The rate given above can be treated as a mean over the redshift range of our sample, 
the mean redshift z = 0.20. We may convert this rate expressed in SNu to the v olumetric rate 



by mu ltiplying the luminosity density. With reference to the luminosity density of iBlanton et al 



(|2003l ). we find the volumetric rate of SNe la, 

rv = (3.63 ± 0.43) x 10"^ Mpc"^ yr'^ (5) 

For g and i passbands the matching of the two LF's yield = 0.278 it 0.036 SNu(g) and = 
0.226 ± 0.034 SNu(i) taking Mg{Q) = 5.07 and Mi{Q) = 4.52 mag, respectively. These SN rates 
yield the volumetric rate ry = (3.86 it 0.50) x 10~^ Mpc~'^ yr~^ for the g passband and ry = 
(4.23 lb 0.64) X 10~^ Mpc~^ yr"^ for the i passband. The rates obtained from the three passbands 
agree with each other, within a 1 a error, showing the internal consistency. These SN rates agree 
with rv = (3.70 ± 0.41 ± 0.34) x 10"^ Mpc'^ yi'^ obtained directly by summing up the LF of SNe 
la in Figure [Tl^a) . Here the second error represents that arising from the completeness function 
given in eq. ([T]), which depends on the SN rate normalisation. If we take the -B-band luminosity 
density Cb = 1-44 x 10*^ Lq Mpc~^ at z = 0.2, obtained by interpolating across five colours, we 
obtain the supernova rate per i?-band luminosity 

TL = 0.257 ±0.028 ±0.024 SNu(B). (6) 

This is compared with earlier measurements at low redshift compiled in Table [3j 

Figure [l2] shows that the SN rate per galaxy is proportional to luminosity of host galaxies, as 
one expects from the comparisons of the luminosity functions. The slope of the solid line is fixed to 
the rate at 0.227 SNu(r) of Figures [9l The data point at the bottom left indicated with the open 
circle is for the 'hostless SNe' where the horizontal error bars indicate only the upper limit on host 
galaxy's luminosity. We conclude that our hostless SNe la are consistent with them being occurred 



The normalisation of the luminosity functions of Blanton et al. 1 20031 ) differs significantly from Blanton et al 



(|200lh apart from the use of the different passbands defined at z = 0.1: the luminosity density of the latter 
20-40% larger than the former. The dif ference in the s hape of the luminosity function, however, is modest and we 



adopt Schechter function parameters of IBlanton et al.l (|200l|) for the luminosity function in the g,r,i passbands at 



z — hy shifting only the normalisation. For the g, r, and i passbands, we interpolated across the five colours at 
2 = 0.1. The luminosity densities at z = 0.2 for the g, r, and i passbands adopted are £g = 1.39 x lO*L0/(Mpc)^ 
Cr = 1.60 X 10®LQ/(Mpc)^ £i = 1.87 x lO*L0/(Mpc)=* 
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in low luminosity galaxies beyond our detection limit but at the rate proportional to luminosity of 
host galaxies. 

We present in Figure [13] the radial distance distribution of SNe la measured from the centre 
of host galaxies measured in the r-passband in physical distance units, where we do not correct 
for the projection effect caused by inclination of host galaxies. The relative distance between SN 
and the centre of h ost galaxy is me asured within the accuracy of 0.1 arcsec, which corresponds to 
0.4 kpc at z = 0.3 (jPier et al.ll2003l ). The ordinate denotes numbers of the SN la per unit volume 
calculated by the formula similar to eq. ([3]). This radial distribution is well represented by the de 
Vaucouleurs profile with the half light radius of = 5.7 kpc as drawn by the thin solid curve in 
Figure [T3l We also draw the exponential profile (dotted curve) with the half light radius of = 3.6 
kpc, which falls off faster than the distribution of SNe la at an impact parameter larger than 10 
kpc. The of the best fit models, = 8.2/11 for the de Vaucouleurs profile and = 11.7/11 
for the exponential profile, where fitting was done within 13 kpc using 13 data points, differ only 
a little but the difference is more apparent in the tail. When one uses the addition of the de 
Vaucouleurs and the exponential profiles, the best fit model shows bulge-to-disk luminosity ratio 
of 0.70, re(deV) = 9.1 kpc and re(exp) = 2.4 kpc with = 5.0/9. This fitted profile is drawn by 
thin dashed curve. 



Our result does not agree with that of iBartunov et al.l (120071 ) who claimed that SNe la in 
spiral galaxies are in a lower rate in the central part compared to SNe la in elliptical galaxies based 
on the ir supernova catalogu e obtained by a compilation of SNe in the literature. Their earlier 
paper (IBartunov et al.lll992l ) claims that the radial dependence of surface density of SNe la can be 
expressed by the exponential profile. Our result does not agree with their conclusion, either. 

The thick solid curve in Figure [13] shows the empirical mean r-passband light profile of field 
galaxies, which is constructed from the aperture flux of galaxies at z = 0.025 — 0.030 in Catalog 
Archive Server of SDSS Data Release 6 using values of profMean. The innermost bin may be 
somewhat affected by finite size seeing. The normalisation is set to the same as Figure [9] while 
taking into acco unt the lumino s ity ev olution of galaxies from z = 0.025 — 0.030 to z = 0.2 using the 
Q parameter of iBlanton et al.l (j2003l ) while fixing the shape of the light profile. The light profile 
of galaxies is consistent with the exponential disc plus de Vaucouleurs spheroid model with the 
bulge-to-disk luminosity ratio of 0.56. We note that the radial distribution of SNe la, at least 
for the bulk of SNe, is consistent with the light distribution (x^ = 11.2/13), except for that at a 
large distance beyond 10 kpc where we see some excess occurrence of SNe. It is interesting to see 
that some SNe la occur at a large distance where the galaxy contributes little light: at > lOkpc 
we expect from the global rate, 1.31 x 10~^ Mpc~^yr~"^ SNe la with the light distribution of the 
de Vaucouleurs profile and 0.15 x 10"^ Mpc~'^yr~^ SNe la with the exponential profile, which are 
compared with observed 0.55 x 10^^ Mpc^'^yr^^ SNe la. The radial distribution of SNe la also 
supports the proposition that the occurrence of SNe la is primarily proportional to the luminosity 
of host galaxies. 
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We study the dependence of the SN rate on colour of host galaxies, to examine whether the 
occurrence of SNe la would correlate with the star formation activity. In Figure [T3] the histogram 
represents the SN host galaxies and the solid curve shows the luminosity weighted colour func tion 
of field galaxies which is calculated using bivariate function (l){Mr,g — r) (jBlanton et al.ll200ll ). as 



^Ug - r)d{g - r) 



'14.75 



dMr(p{Mr,g- r) x lO°-<*^©.'-^^'-)(i(5 - r), 



(7) 



24.25 



where (p{Mr,g- r) is evaluated for -24.25 < Mr < -14.75 and 0.12 < g-r < 0.88 on the 20 x 20 
grid. The abscissa is rest-frame colour after the K-correction. The normalisation is taken to be 
the same as that in Figure El Morphologic al types of galaxies are indicated with the corresponding 
colours according to iFukugita et al.l ()2007l ) with the error bar representing the dispersion of colours 
in the morphologically classified sample. This figure shows that the colour distribution of SNe la 
host galaxies traces well that of field galaxies, and we do not see any particular excess of SNe la 
host galaxies for bluer, late type or irregular galaxies, beyond one sigma level. In this comparison 
we do not take into account the evolution of galaxy colour to z = 0.2. 

In particular we may be interested in the difference of the SN rate between the elliptical 
galaxies and other galaxies. If we select galaxies with colour 0.73 < g — r < 0.81 corresponding 
to elliptical galaxies (with contaminations from SO and some Sa galaxies) we obtain the SN rate 
0.194ib0.062 SNu(r) which is compared with 0.255ib0.035 for galaxies with g — r < 0.73 (Sa galaxies 
or later). We observe a 31±35% enhancement in the SN rate relative to galaxy luminosity in late 
type galaxies compared to that in elliptical galaxies, though the effect is only at one sigma. We 
do not see a particular enhancement of SN rate per luminosity, however, among late type galaxi es, 
where the star formation rate increases towards later type morphologies (jNakamura et al.ll2004l ). 



u — g colour is more sensitive to the star formation activity than g — r, but u-band brightness 
is too faint to us for most host galaxies: SN host galaxies brighter than 21 mag are less than 40% 
for our magnitude limit supernova sample, and we are not able to draw a meaningful conclusion 
from the u — g colour distribution. The SN host galaxies are also too faint for SDSS spectroscopy 



(jStrauss et al.ll2002l ). and we cannot estimate the star formation rate directly, unless we resort to 
population synthesis colours, which are not well constrained without u colour. Therefore, we do 
not find clear evidence that points to the correlation of the SN rate with the star formation activity. 



Mannucci et al.l (|2005l ) and I Sullivan et al.l (|2006l ) give models of SN la rates with the two com- 



ponents of progenitors, one explosions of old binary stars with a large delay t ime from the for r nation 
and the other explosions with a short delay time after the system formed. iMannucci et al.l ()2005l ) 
give SNe la rate that depends on the morphological type of host galaxies in units of SNu(K). 
We may use their rate to estimate the number of SNe la for each morpholo gical type using 
morph ological-type dependent luminosity density of field galaxies, as given by 



Nakamura et al 



(j2003l ) with the aid of the conversion from the r to the X-passband (iNagamine et al.ll2006l ). This 
model gives ry = 3.0 x 10^^ Mpc~^ yr^^, where 31% are the prompt component, which is consis- 
tent with the difference in SN la rate relative to galaxy luminosity between elliptical galaxies and 
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spiral galaxies. The total rate also agrees with our rate in equation ([5]) 
presented in Table [H 



The detailed numbers are 



Sullivan et alj (|2006l ) modelled SN la rate as a function of stellar m a,sses and mean SFRs o f 



host galaxies. We use the morphology dependent luminosity density of iNakamura et al.l (l2003l ) 
and the morphologically dependent star fo rmation rate cal c ulated from the star formation rate 
of the bugle and disk components given in iNagamine et al.l (|2006l ). The predicted SN la rate is 
ry = 2.6 X 10"^ Mpc^^ yr^^j of which 13% are prompt: see also Table [H The predicted numbers 
of SNe la are plotted in Figured!] as dotted (Mannucci et al.'s model) and dashed (Sullivan et al.'s 
model) curves taking t he colour distribution of each morphological type as a Gaussian with mean 
and dispersion given bv lFukugita et al.l ( 20071 ). The for the two models are 15.1 and 23.6 for 11 
degrees of freedom using data for 0.1 < g — r < 1.0. 

The currently available two component models are marginal to represent our data, but are 
consistent with the observed number of SNe la, as shown by the histogram in Figure [T4l allowing for 
large uncertainties. The prompt component is rather minor and not very manifest in the integrated 
rates for low z SNe, as in our analysis. The model prediction for Im galaxies is low, which is due 
to the low luminosity density of Im galaxies in the SDSS morphologically classified sample. In 
Nakamura et al.l (120031 ) Im galaxies contribute only ~ 1% of the total luminosity density, which 
may be due to the sel ection effect that disfa vours Im galaxies. This contrasts to the indication from 
the colour function of lBlanton et al.l (|200ll ). which suggests that galaxies bluer than {g — r) < 0.42 
contribute by ~ 15% of total luminosity. If the luminosity density of Im galaxies is this large, 
however, the number of SNe la in Im galaxies would be 10 times more that would largely overshoot 
our observed numbers of SNe. To identify the prompt component and test the two component 
hypothesis, we need a sample where the prompt component is dominant, or the sample with which 
both star formation rate and SN rate can more accurately be measured. 

Our sample allows us to study if the properties of SNe la would depend on the property of host 
galaxies. It has occasionally been claimed that bright SNe la appear more often in late type galaxies , 



while SNe la in ea r ly type galaxies are subluminous (jFilippenko 



198S 



della Valle Sz Panagia 



1992 



Reindl et al.l l2005l : ISullivan et al.l l2006l ) . This has been counted as one of the reasons that SN la 
calibration of nearby galaxy's distance led to a smaller value of the Hubble constant in the past. 
Figure [TSTa) plots the maximum brightness of SNe la in our sample as a function oi g — r colour of 
host galaxies: mean brightness of SNe la does not change across colours of galaxies from E to Im. In 
particular the plot does not indicate any evidence that SNe la in bluer galaxies are systematically 
brighter or those in red galaxies are fainter. Both mean and dispersion of luminosity are nearly 
constant with respect to galaxy colour as seen in the figure. 

Figure [TSlfb) shows the plot of the xi parameter, which represents the decline rate of SNe la 
brightness, against g — r colour of host galaxies. The upper values of xi does not change with 
respect to colour from E to Im galaxies, but the lower value shows the trend that it decreases from 
—0.5 (or Ami5 ~ 1.2) for Im galaxies to —3 (Am-15 ~ 1.7) for E and SO galaxies. It is noted that 



- 13 - 



large negative values of xi are seen only in early-type galaxies. This causes some dependence of xi 
on colour of host galaxies as seen in the figure. The slope, however, is small enough to cause any 
effect on brightness of SNe la through the correction of 0.12x1. 

Similarly Figure [16] shows maximum brightness of SNe la on luminosity of host galaxies, 
showing that the former does not depend on luminosity of the host. Note that we expect that 



avera ge metallicity of the host galaxy changes by 0.59 dex in this luminosity range ([Tremonti et al 



2004 ) . The xi parameter changes slightly as luminosity of host galaxies changes. 



We also examine the dependence of the c parameter and the xi parameter on the distance 
from the centre of galaxies as shown in Figure [TTl The colour excess c = E{B — V) stays nearly 
at ~0.08 with the dispersion of 0.07 and does not show a systematic change from 1 kpc to 15 kpc. 
In particular, we observe that the colour excess does not decrease as we go farther away from the 
centre of galaxies, where stars, and hence dust, are expected to decrease. The trend we see here 
does not change if we limit the SNe to z < 0.2 with which SNe la that would receive reasonably 
large extinction are included in the sample. This finding suggests us to interpret that the observed 
reddening is mostly associated with individual supernovae, either extinction from supernova itself 
and/or circum-supernova dust or intrinsic colour variation of the supernova, rather than interstellar 
dust in host galaxies. This would also give us an upper limit on the model as to the amount of 
dust produced around SNe la. We observe only 7 SNe la which shows colour excess more than 
E(B — V) > 0.3 mag among our 137 SN la sample. In addition, we do not observe a systematic 
dependence of the xi parameter on the distance from the centre of galaxies; see Figure fTTlf b). 



6. Conclusion 

The SNe la sample acquired in the SDSS II, containing 137 low redshift SNe from z = 0.05 to 
0.3, indicates that the occurrence of SNe la is primarily proportional to luminosity of galaxies. The 
LF of SN la host galaxies matches very well with that of field galaxies multiplied by luminosity, and 
colour of SN la host galaxies does not differ from that of field galaxies. Our low redshift sample 
does not indicate an active signature that the occurrence of SNe la follows star formation activity, 
except that possible enhancement in the SN la rate is noted in late type galaxies (< 31 ± 35%) 
compared with the rate in elliptical galaxies. Our low redshift sample is compatible with the two 
component model where the effect of the prompt component is modest (10—30% of the total rate), 
as in the current models. We are not able to differentiate SN rates among late type galaxies. Our 
sample contains 8 SNe la, whose host galaxies were not identified. It is shown, however, that they 
are consistent with them occurred in low luminous galaxies beyond the survey limit for galaxies. 
Luminosity of SNe la does not appear to depend upon luminosity or colour of host galaxies. 

The luminosity function of SNe la is Gaussian with Mb = —19.42 and a = 0.24 mag (FWHM is 
a factor 1.4 in luminosity), if the colour variation is interpreted as reddening obeying the extinction 
law of the Milky Way with the standard value Ry = 3.1. This Gaussian distribution is further 
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tightened and a = 0.14 mag if the extinction to reddening ratio Ry is reduced to ~ 2 and if the 
correlation is taken into account between maximum brightness and the decHne rate. Reddening 
of SN la, inferred from the colour variation, is 0.1 mag in Ay, which is lower approximately by a 
factor of 5 than that for star forming field galaxies measured from Ha and H/? emission. 

We also showed that the distribution of the SN la occurrence, that extends to a few tens of 
kpc, agrees with the general light distribution of galaxies that have bulges with the de Vaucouleurs 
type profiles. The variation of colours of SN la is constant and does not depend on the distance 
from the centre of galaxies unlike what is expected for the supernova that would happen in the 
galactic disc. The colour variation of SNe la may likely be ascribed to intrinsic colour variation or 
immediate neighbourhood of supernovae rather than extinction in host galaxies. The total column 
density of dust should not give extinction more than 0.2 mag in the V passband. 
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A. Verification of the analysis method 

In order to examine the validity of our analysis, we have made simulations for our observation 
creating a set of light curves using SALT2 model in combination with SNANA. The distribution of 
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xi parameter is assumed to be a Gaussian with the mean and the dispersion 0.90. The distribution 
of c is a Gaussian with the mean 0.063 and the dispersion 0.11. We assume the relation involving 
absolute magnitude M, xi and c, 

M = Mo-axi+ (3c (Al) 

with a = 0.10 and /? = 2.45. The intrinsic dispersion of 0.11 mag was applied for every passbands 
and every epochs. The epochs of observation is fixed to actual date and the signal to noise ratio 
of each observing point was calculated by using real observational condition of our observation 
including seeing and sky brightness. The detection efficiency searching for variable objects is also 
included in our simulation, and the SN la rate of rsNeia = 2.2 x 10^^ (1 + zY'^ (Mpc^^yr^^) is 
assumed. The cut based on the quality of light curves as described in the text is also applied to the 
simulation, and the spectroscopic incompleteness of eq. ([T]) is applied with Zc = 0.15 and Az = 0.3. 
The simulation generated 210 SNe light curves. We make 50 sets of simulated light curves and 
applied the method the same as that to the observed dataset. The incompleteness correction e{z) 
in eq. ([3]) is calculated separately for each simulation set. 

Figure [18] shows the comparison of the input and the output SN LF after correcting for the 
effect of c parameter. Solid curve represents the expected LF taking into account the error of 
the determination of and c from light curve fits. Histograms and error bars are the mean 
and dispersion in each magnitude bin calculated from the 50 simulated light curve sets. A good 
agreement is seen between solid curves and histograms. 

We also make a simulation incorporating the completeness as a function of apparent peak 
r-passband brightness of SNe like e(m) = 1.0 — (rmax — 19.75)/3.15 to see the effect when the 
completeness is not a function of redshift but of apparent peak magnitude. The completeness e{z) 
in eq. ([3]) is estimated as a function of redshift. The result is shown in Fig. [T9j We see again 
a good agreement between the input and output SNe. This is due to the fact that redshift and 
apparent peak magnitude are well correlated. 

Good agreements are seen between input and output SN LF's for both cases where incomplete- 
ness is a function of redshift and of apparent peak r-passband magnitude. The analysis adopted 
in this paper does not suffer from particular systematic effects while dealing with apparently a 
dataset with sample incompleteness. The volumetric SN rate is recovered to be (2.95 ±0.42) x 10~^ 
and (3.04 ± 0.51) x 10~^ depending on whether incompleteness is taken to be a function of red- 
shift or apparent magnitude, respectively. These numbers are compared with the input value of 
3.18 X 10~^ (Mpc~^ yr^^). The volumetric SN rate is properly recovered within the error. 

We also examine the distribution of the c and xi parameters of SALT2 light curve fit in 
Fig. [20l SNe with high c value or low xi value will be faint and such SNe may suffer from 
incompleteness. The analysis is similar to that of the LF. In eq. [3l we replace absolute magnitude 
M with the colour parameter c or the shape parameter xi. We do not see any systematic effect 
for the recovered distribution. We change the assumed distributions of the c and xi parameters 
and the relation among M, xi and c in a way they are consistent with the observed SNe, and 
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calculate the completeness of the observed sample. This procedure is iterated till the convergence 
(one or two iterations are sufficient). The final values are xi = i: 0.96, c = 0.039 it 0.093 and 
M = Mq — 0.12xi + 2.52c, where the last relation agrees with the one obtained by minimising the 
width of the LF of SNe la as seen in Fig. E^d). With these procedures the redshift distribution of 
the observed SNe is reproduced as shown in Figure [2TJ 

We conclude that the method we used in this paper works well even for the data with incom- 
pleteness after the proper account taken for incompleteness. We do not expect any specific biases 
in the analysis. 
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Table 1. Sample selection summary 



Criteria Number 



total 314 

spectroscopic/confm'd 130 

spectroscopic/probable 16 

light curve selection 168 

secure light curve 222 

good sky area date 207 

riim < 21.5 137 

no host galaxies 8 



Table 2. Hostless SNe la 



SDSS ID 


lAU Name 


Redshift 


Luminosity Limit (L*) 


Dist. of 10"(kpc) 


2943 


2005go 


0.2659 


0.055 


40.7 


5994 


2005ht 


0.1885 


0.025 


31.2 


6780 


2005iz 


0.2046 


0.029 


33.1 


6924 


2005ja 


0.3286 


0.089 


47.2 


6933 


2005jc 


0.2137 


0.033 


34.5 


7475 


2005jn 


0.3188 


0.085 


46.7 


7335 


2005kn 


0.1975 


0.028 


32.6 


3565 




0.2885 


0.067 


43.4 


8030^ 


2005jv 


0.4226 


0.161 


55.5 



^Outside the magnitude limited sample 
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Table 3. SNe la Rate Measurements 



Reference Redshift Rate (SNu) 



Cappellaro et al. (1999) 





0.18 ±0.05 


Dilday et al. (2008)^ 


0.09 


246+°-°^^ 

U.Z4D_Q QgQ 


Madgwick et al. (2003) 


0.098 


0.196 ±0.098 


Blanc et al. (2004) 


0.13 


n 1 oc;+0.044+0.028 
U.lzo_Q Q34_Q 


Hardin et al. (2000) 


0.14 


n 99+0.17+0.06 
"•^"^-0.10-0.03 


This work 


0.20 


0.257 ±0.028 ±0.024 


Botticella et al. (2008) 


0.3 


r, 99+0.10+0.16 
^•^^-0.08-0.14 



^B-hand luminosity density of js = 1-19 x 10*^ Lq Mpc 
was used to convert from the value per comoving volume unit. 



Table 4. Predicted SNe la rate in the two component model 













Mannucci et al. model 






Sullivan et al. model 




morph. type 


g-r 


jr 




Na (delayed) 


Nb (prompt) 


Na + Nb 


Na (delayed) 


Nb (prompt) 


Na + Nb 






loSLr.(0)Mpc-3 




lO-^Mpc-^yr-i 


lO-^Mpc-^yr-i 


lO-^Mpc-V"^ 


lO-^Mpc-V"^ 


lO-^Mpc-^yr-i 


IQ-^Mpc-^yr-i 


E,SO 


0.75 


0.62 


3.32 


0.720 


0.000 


0.720 


0.882 


0.052 


0.935 


SO/a-Sb 


0.64 


1.00 


2.77 


0.970 


0.304 


1.274 


1.052 


0.189 


1.240 


Sbc-Sd 


0.51 


0.37 


2.52 


0.326 


0.495 


0.821 


0.326 


0.088 


0.414 


Im 


0.36 


0.02 


2.15 


0.015 


0.127 


0.142 


0.013 


0.006 


0.019 


Total 




2.01 




2.031 


0.926 


2.957 


2.273 


0.335 


2.608 



I 

tN3 
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Redshift 



Fig. 1. — Redshift distribution of our samples. Tlic dotted liistogram represents the basic sample 
comprising 207 SNe la. The solid histogram is the magnitude limited sample, our final sample, 
comprising 137 SNe la with rmax < 21.5 mag used to derive the LF. Shaded region shows the 
expected distribution when SNe la are distributed as oc (1 + z)^-^ with the selection criteria we 
adopted taken into account. Shading stands for 20% uncertainty in the normalization. Dashed 
curve is the expected distribution when SNe la rate shows no evolution in z. 
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Fig. 2. — Apparent maximum brightness of 207 SNe la in the r-band as a function of redshift. 
Horizontal dotted Hne denotes the magnitude hmit of r^j^ = 21.5 mag. Brightness is not corrected 
for extinction or colour excess. 
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Fig. 3. — Solid line shows the ratio of the number of SNe la in our basic sample to the simulation 
as a function of redshift. Dotted line shows the completeness function eq.([I]) used in this paper. 
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Fig. 4. — Distributions of apparent r-band Pctrosian magnitude of host galaxies of our basic SN 
la sample (dotted histogram) and of the magnitude limited sample (solid histogram) . There are 9 
(or 8) SNe la whose host galaxies are not identified in the two samples, which are indicated in the 
leftmost bin with shades. 



-26- 



1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



0.8 - 



2 0.6 

X 

m 
I 

O 

a 



0.4 - 



0.2 



A= 0.048±0.012 
G= 0.074±0.010 





-0.5 



J L 



I I I I 



0.5 
E(B-V) = c 



Fig. 5. Distribution of colour excess E(B — V) = c calculated from the magnitude limited sample. 
The curve is the Gaussian with the function p{E{B — V)) oc exp{—E{B — V)/ A) for E(B — V)>0 
and p{E{B — V)) = 0.0 for E(B — V)<0 convolved with Gaussian (see text). 
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Fig. 6. — Distribution of y-band extinction Ay for SNe la, estimated assuming that their colour 
excess is due to dust reddening and Ry = 3.1 (solid histogram). T he dotted histogram is Ay 
of field galaxies inferred from the Balmer emission line ratio Ha/H/3 ([Nakamura et al.l |2004| ) for 
comparison. 
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Fig. 7. — Luminosity function of SNe la in the 5-band. B-band brightness is corrected for the 
colour variation as Mb — /3c (/3 = Ry + 1) with (a) /? = 4.1 and (b) /? = 2.92 which minimises 
the width of fitted Gaussian. The parameters of the Gaussian are shown in each panel. In panels 
(c) and (d) a correction for light curve shape X\ is also taken into account with the parameter to 
minimise the Gaussian widths. 
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Fig. 8. — Correlations between brightness of SNe and the colour excess parameter c (the top panel) 
and the xi parameter (the bottom panel). 
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Fig. 9. — Luminosity function of SNe la host galaxies in the r-passband. Lower panel shows the 
number of contributing ga laxies in each bin. T he solid curve is the luminosity weighted luminosity 
function of field galaxies ([Blanton et al.ll200ll ) normalized to fit the histogram. The scale factor 
that gives the supernova rate in SNu is given at the top right of the figure. Dotted lines show 
the range corresponding to the fitting error. Contribution from hostless SNe are indicated in the 
rightmost bin with shades. 
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Fig. 10. — Same as figure [9] but in the ^-passband. 
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Fig. 11. — Same as figure [9] but in the i-passband. 
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Fig. 12. — SN la rate per galaxy as a function of r-band luminosity of galaxies. The point denoted 
by the open symbol stands for the hostless SNe la and the horizontal bar shows the upper limit of 
luminosity of galaxies. 
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Fig. 13. — Radial distribution of SNe la from the centres of their host galaxies measured in physical 
distance scale. The thin solid, dotted and dashed lines are best fit models with de Vaucouleurs, 
exponential profiles and sum of these two profiles, respectively. The thick solid line is the mean 
r-passband galaxy light profile averaged over field galaxies. 
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Fig. 14. — g — r colour function of SNe la host galaxies with the lower panel indicating the number 
of contributing galaxies in each bin. Solid curve represents the luminosity weighted galaxy colour 
function normalized in the way same a s in Figure [9l M o rphol ogical types of galaxies are indicated at 
the corresponding colours according to lFukugita et al.l (120071 ). Horizontal bars indicate the variance 
within the morphological ty pes. Dotted and dashe d lin es are numbers expec ted for SNe la based 
on the SNe la rate model of lMannucci et al.l (j2005l ) and ISullivan et al.l (j200a ). respectively. 
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Fig. 15. — (a) Absolute brightness of SNe la corrected for the colour excess using /? = 4.1 as a 

function of {g — r) colour of host galaxies. The mean and dispersion of the points arc also plotted, 
(b) Shape parameter xi that depends on colour (g — r) of host galaxies. Dotted lines show least 
square fits to the data, Mb - 4.1c = -19.32 + 0.01(5 - r) for (a), and xi = 0.26 - 0.74(5 - r) for 
(b), respectively. 
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Fig. 16. — (a) Absolute brightness of SNe la corrected for the colour excess using (5 = 4.1 and (b) 
shape parameter x\ as a function of r-band brightness of host galaxies. Dotted lines show least 
square fits to the data, Mb - 4.1c = -19.78 - 0.034Mr for (a), and xi = 3.50 + O.ISM^ for (b), 
respectively. 
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Fig. 17. — Colour excess and the xi parameters plotted against the physical distance of SNc la 
from the centre of host galaxies. Dotted lines indicate the zero level. In upper panel the mean and 
dispersion in respective bins are also plotted. 
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Fig. 18. Comparison of the input SN LF (solid curve) and output SN LF (histogram) from 

simulated datasct assuming the completeness as a function of redshift. Histogram passes through 
the mean and error bars show the dispersion in each bin calculated from 50 simulation dataset. 
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Fig. 19. — Same as Fig. [18] but completeness is assumed to be a function of apparent peak r-band 
magnitude. 
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Fig. 20. — Same as Fig. [18] but for the distribution of the E(B — V) = c parameter (left) and the 
xi parameter (right). 
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Fig. 21. — Simulated redshift distributions before (solid) and after (dotted) spectroscopic incom- 
pleteness is applied. Points with error bars represent observed SNe sample. 



